SHORT-PULSE  RADAR  VIA  ELECTROMAGNETIC  WAVELETS 


Proceedings  of  the  Third  Ultra- Wideband,  Short-Pulse  Electromagnetics  Conference 

Albuquerque,  NM,  May  27-31,  1996 
To  be  pubhshed  by  Plenum  Press 


Gerald  Kaiser 

Department  of  Mathematical  Sciences 
University  of  Massachusetts  Lowell 
Lowell,  MA  01854,  USA 
Email:  gkaiser@cs.uml.edu 


1 9970527  023 


INTRODUCTION 

The  new  theory  of  physical  wavelets  makes  it  possible  to  perform  radar  and  sonar 
analysis  directly  in  the  space-time  domain,  based  on  fundamental  principles  underlying 
the  emission,  reflection,  and  reception  of  electromagnetic  and  acoustic  waves^’^.  Being 
independent  of  the  Fomrier  transform  and  even  of  the  usual  (affine)  wavelet  transform, 
this  formalism  is  therefore  equally  well  suited  for  ultrawideband  or  short-pulse  radar  as 
for  narrowband  or  continuous-wave  radar.  However,  Fourier  analysis  does  have  a  natural 
place  in  this  theory  and  can  be  used  easily  when  spectral  questions  are  of  interest. 

A  transmitting  antenna  following  an  arbitrary  (possibly  accelerating  or  nonlinear) 
space-time  trajectory  a{t)  emits  a  physical  (acoustic  or  electromagnetic)  wavelet  which  is 
propagated  in  space  by  the  appropriate  Green  function.  This  defines  an  emission  operator 
Ea  which,  acting  on  any  time  signal  ip{t),  gives  the  emitted  space-time  wave  {Ea^){T,t). 
The  reception  operator  Ra  is  di’al  to  Ea  ,  measuring  any  incident  space-time  wave  F(r,  t) 
along  the  given  antenna  trajectory  a(t)  to  produce  the  received  time  signal  (RaE)(t). 
Reflection  is  modeled  as  reception  followed  by  re-emission,  i.e.,  by  the  operator  EaRa 
transforming  any  incident  space-time  wave  to  the  reflected  space-time  wave.  Let  the  re¬ 
ceiving  antenna  follow  another  arbitrary  space-time  trajectory  7(t)  (possibly  different  from 
the  trajectory  a(t)  of  the  transmitter),  let  the  target  follow  a  third  arbitrary  space-time 
trajectory  /3t(^)>  let  the  transmitted  and  received  signals  be  %p{t)  and  x{i)i  respectively. 
The  objective  is  to  estimate  the  target  trajectory  Prif)  from  a  knowledge  of  Q(t),  'y{t),ip{t) 
and  x(i)-  This  is  achieved  by  maximizing  the  modulus  of  the  the  normalized  ambiguity 
functional  Xn{P),  obtained  by  matching  the  actual  return  x(f)  with  the  computed  return 
due  to  a  trial  trajectory  P{t).  When  the  radar  is  monostatic  and  the  target  is  assumed  to 
move  uniformly  in  the  radial  direction,  then  Xn{P)  reduces  to  the  usual  wideband  ambi¬ 
guity  function,  which  is  just  the  ordinciry  time-scale  (wavelet)  transform  of  x(i)  with  i/^(t) 
as  the  basic  wavelet.  In  the  narrowband  approximation,  it  reduces  further  to  the  usual 
time-frequency  ambiguity  function,  which  is  a  windowed  Fourier  transform  of  the  video 
signal  of  the  return.  This  shows  that  our  “physical  wavelet  analysis”  is  a  generalization 
of  the  usual  (“mathematical”)  wavelet  analysis,  which  is  in  turn  a  generalization  of  time- 
frequency  analysis.  In  particular,  our  analysis  applies  equally  well  to  bistatic  wideband 
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radar,  where  the  Doppler  effect  can  no  longer  be  represented  simply  by  scaling  and  hence 
the  usual  affine  wavelet  analysis  breaks  down. 

In  Reference  2,  the  transmitting  and  receiving  antennas  and  the  target  were  all  as¬ 
sumed  to  be  points,  so  that  the  trajectories  a{t),l3{t),  and  7(t)  fully  describe  their  motions. 
Consequenctly,  the  emitted  and  reflected  waves  are  omnidirectional  and  hence  are  not  very 
useful  in  practice.  In  this  paper  we  generalize  the  above  model  to  include  extended  an¬ 
tennas  and  targets,  follwing  an  idea  introduced  by  Heyman  and  Felsen^.  The  associated 
radar  beams  are  much  more  useful  since  they  have  a  measure  of  directivity. 


EXTENDED  PHYSICAL  WAVELETS  AND 
AMBIGUITY  FUNCTIONALS 


Suppose  we  are  given  an  antenna  located  at  the  space  point  x,  emitting  the  response 
to  an  impulse  at  time  t.  Ignoring  polarization  for  simplicity,  the  resulting  wave  at  the 
observation  point  x'  at  time  t'  can  be  represented  by  a  solution  of  the  scalar  wave  equation, 
which  we  write  as 

K{x!,t'\x,t), 

with  a  source  distribution  appropriate  to  the  antenna.  If  the  antenna  is  very  small  (essen¬ 
tially  a  point)  and  omnidirectional,  then  A(x',  t'  \  x,  t)  is  well  approximated  by  the  retarded 

Green  function  ,  ,  ^ 

.  ,  ,  s  t  -  |x' —  x|/c) 

- 4^1^,  _  ^1 -  . 


where  c  is  the  speed  of  light.  Following  Heyman  and  Felsen^,  a  simple  model  can  be 
formulated  for  an  extended  antenna  by  allowing  complex  antenna  space-time  coordinates 
X  z  =  X  +  fy  and  t  -^u-t  +  is  and  taking  K{x',  t'  \  z,  u)  to  be  an  analytic  extension 
of  the  above  retarded  Green  function.*  Heyman  and  Felsen  showed  that  for  |y|  <  cs, 
K{x' ,t' \z,u) .  can  be  interpreted  as  a  pulsed  beam  Geld  emitted  by  a  circular  disk  of 
radius  |y|  in  the  direction  of  y.  Thus  y  is  a  convenient  “handle”  by  which  the  radius  and 
orientation  of  the  antenna  can  be  controlled,  without  having  to  construct  a  messy  model 
for  the  antenna  involving  a  continuous  distribution  of  point  sources.  More  complicated 
extended  antennas  and  arrays  can  be  modeled  by  a  distribution  (continuous  or  discrete) 
of  such  complex  source  points. 

To  keep  the  notation  uncluttered,  we  combine  the  space  and  time  coordinates  into  a 
single  symbol: 

x'  =  (x^^')eR^  z  =  {z,u)  =  x  +  iyeC^, 


where 

x  =  {x,t)  and  y  =  {y,s). 


*  Note  that  our  convention  differs  from  that  of  Reference  3  in  that  our  positive- 
frequency  time-harmonic  waves  vary  as  rather  than  e  .  For  this  reason,  the  analytic 
continuation  of  the  retarded  Green  function  (using  the  analytic  signal  of  the  delta  func¬ 
tion)  is  to  the  upper-half  time  plane  (s  >  0)  rather  than  the  lower-half  time  plane.  This  is 
consistent  with  the  convention  used  in  Reference  1. 


2 


r 


The  condition  |y|  <  cs  means  that  the  imaginary  space-time  four-vector  y  belongs  to  the 
future  cone,  so  that  z  actually  belongs  to  the  complex  future  tube^’'^ 

z  £  T+  =  {x  +  iy  £  :  X  ^  and  y  =  (y, s)  with  |y [  <  cs}, 

which  is  a  four-dimensional  generalization  of  the  upper-half  complex  time  plane. 

Suppose  now  that  om  antenna  executes  an  arbitrary  motion,  including  possible  rota¬ 
tions  and  accelerations.  Using  the  complex  source  coordinates,  this  can  be  parameterized 
as 

z  =  a{t)  =  x{t)  +  iy{t),  where  x{t)  =  (x(t),  t)  G  R"*  and  y{t)  =  (y(t),  s).  (1) 

It  is  reasonable  (but  mathematically  unnecessary)  to  assume  that  the  radius  of  the  antenna 
remains  constant  dming  the  motion,  so  that  |y(t)l  =  ly(0)|  =  R  <  cs,  although  the 
direction  of  y(t)  may  vary  to  allow  tracking,  scanning,  etc.  While  executing  this  motion, 
the  antenna  is  fed  an  input  time  signal  ip{t).  Then  the  output  beam  is 


dt  K{x'  I  a{t))ip{t). 


For  reasons  explained  in  Reference  2,  we  call  the  extended  physical  wavelet  gener¬ 

ated  by  ip(t)  along  the  antenna  motion  a(t).  Given  a(t),  we  define  the  emission  operator 
Eq  as  the  operator  transforming  the  time  signal  '^(t)  to  the  space-time  wave  i.e., 

/OO 

dtK{x'\a{t))^{t).  (2) 

-OO 


Thus  Eq  takes  a  function  of  one  variable  (the  input  signal)  to  a  function  of  four  variables 
(the  output  beam).  On  the  other  hand,  if  the  antenna  is  used  as  a  receiver,  it  converts 
space-time  waves  into  time  signals.  Again,  assume  that  the  complex  antenna  motion  a{t) 
is  given  as  in  (1).  Then  the  simplest  model  for  the  received  signal  due  to  an  incident  wave 
F(x')  is 

{RQF){t)=gQF{a{t)),  (3) 

where  Qq  is  a  “gain  factor.”  Thus  Rq  simply  measures  the  field  along  the  complex  tra¬ 
jectory  a{t).  More  complicated  receivers  can  be  formulated  which  measure  derivatives  of 
F  along  a{t).  (In  the  full  electromagnetic  formalism,  for  example,  Rq  could  measure  the 
induced  current  rather  than  the  field.)  Since  a(t)  is  complex,  the  “evaluation”  of  the  field 
F(x')  at  x'  =  a{t)  must  be  defined  in  (3).  For  this  we  use  the  analytic-signal  transform  of 
F,  which  extends  F  to  complex  space-time^’^: 

r./  sip  dr 

F{x  +  iy)  =  —  /  - :  F{x-\-Ty). 

J-oo  T'  -  i 

When  y  =  (0,  s)  with  s  >  0,  F{x-\-iy)  reduces  to  the  usual  Gabor  analytic  signal  F(x,  t+is) 
corresponding  to  F(x,  t),  with  x  regarded  as  an  external  parameter;  this  function  is  analytic 
in  the  upper- half  complex  time  plane.  It  is  further  shown  in  Reference  1  that  if  F{x')  is  any 
solution  of  the  homogeneous  wave  equation  (or  Klein-Gordon  equation^),  then  F{x  -t-  iy) 
is  analytic  in  the  future  tube  7^  (i.e.,  ly|  <  cs).  The  reception  operator  then  evaluates  the 
analytic-signal  transform  F{x  +  iy)  in  its  region  of  analyticity. 
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With  emission  and  reception  modeled  by  (2)  and  (3),  we  are  almost  ready  to  formulate 
a  general  radar  problem.  The  only  missing  element  is  a  model  for  reflection.  In  the  spirit 
of  regarding  a  scattered  electromagnetic  wave  as  being  emitted  by  the  current  induced 
on  the  scatterer  by  the  incident  wave,  we  propose  the  following  model:  Suppose  we  are 
given  an  oriented  circular  “target”  disk  executing  a  motion  described  by  a  complex  space- 
time  trajectory  a(t)  =  x{t)  +  iy{t)  as  in  (1).  Again,  we  interpret  the  imaginary  position 
vector  y{t)  as  defining  the  radius  and  orientation  of  the  disk.  (To  say  that  the  disk  is 
“oriented”  means  that  its  two  sides  are  not  equivalent;  for  example,  one  side  could  be 
reflective  while  the  other  side  is  not.  Then  every  unit  vector  y  =  y/ly|  corresponds  to  a 
unique  orientation  of  the  disk.  This  is  useful  if,  for  example,  we  approximate  a  complicated 
target  by  patching  together  disks  of  various  sizes  and  orientations,  as  in  Section  10.2  of 
Reference  1;  their  non-reflecting  sides  should  then  be  oriented  towards  the  interior.)  A 
given  space-time  wave  F{x')  will  now  be  assumed  to  be  reflected  from  the  disk  as  follows: 
First  the  disk  acts  as  a  receiver,  then  as  a  tra,nsniitter .  Thus  the  reflected  wave  is 

/CO 

dt  K{x'\cx{t))F{a{t)). 

-OO 

Note  that  in  the  present  context,  the  original  “gain  factor”  Qa  is  re-interpreted  as  a.  re¬ 
flection  coefficient.  When  a  complicated  target  is  patched  together  from  circular  targets 
of  various  radii  and  orientations,  the  reflection  coefficient  becomes  a  function  defined  over 
the  target  surface  as  desired. 

The  ambiguity  functional  formalism  developed  in  Reference  2  generalizes  easily  and 
naturally  to  the  present  setting  of  extended  physical  wavelets.  Given  the  outgoing  time 
signal  V’  the  motions  a,  P ,  and  7  of  the  transmitter,  target,  and  receiver  (all  complex), 
our  model  for  the  time  signal  received  at  7  is 


M^")  =  {RiEpRpEamt") 

r  r  (4) 

= s-yso  /  /  dt'  dt  K(i(t")  I  m)  KW)  I  «(<))  m- 

Of  course,  the  received  signal  depends  functionally  on  all  three  trajectories  cx.,P,^,  as  is 
evident  from  the  right-hand  side  of  (4).  But  to  simplify  the  notation,  we  have  suppressed 
the  dependence  on  the  known  trajectories  a  and  7  and  displayed  only  the  dependence  on 
the  target  trajectory  p.  To  estimate  the  actual  target  trajectory  Pr{t),  we  compute  'ip0{t) 
for  a  trial  trajectory  p{t)  and  match  the  result  with  the  actual  return  x{t)  by  taking  the 
inner  product  of  the  two  time  signals.  We  denote  the  result  by  x{P)t  which  we  call  the 
ambiguity  functional  of  the  return: 


/OO 

dt"  x(t") 

-OO  (5) 

-III  dt"  dt'  dt  x{t")  I  Pit'))  KiPit')  I  ait))  iPit). 


(We  assume  that  V'(t)  and  xit)  are  real;  if  they  are  complex,  then  x(t)  should  be  replaced 
by  its  complex  conjugate  in  (5).)  Assuming  that  Mt)  and  xit)  have  finite  energies 
and  11x11^,  the  Schwarz  inequality  implies  that 


lx(;0)!  =  \ix,i^i3)\  <  11x11  W'fppW 
\xiP)\  =  11x11  Wi’isW  xit)  =  Cippit). 
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Therefore,  to  estimate  the  true  target  trajectory  /?r(0)  maximize  the  normalized 

ambiguity  functional 


Xn(/5)  = 


m 
m\  ■- 


By  (6), 


IXn(/5)  I  <  llxll  and  \xn{P)\  =  11x11  x{t)  =  Cipp{t). 


Equivalently,  we  can  minimize  the  error  functional  defined  by 


e{P)  =  1  - 


\xm 
llxll  UpV 


since  the  Schwarz  inequality  states  that 

0  <  S{P)  <  1  and  £{P)  =  0  x(^)  =  C'ipp{t). 

Thus  |xAf(/?)|  and  E{0)  attain  their  maximum  and  minimum  values,  respectively,  only 
when  the  trial  retmn  is  indistinguishable  from  the  actual  return.  Of  course,  this  does  not 
guarantee  that  the  trial  trajectory  fi{i)  coincides  with  the  actual  target  trajectory  (drit)-, 
since  the  return  does  not,  in  general,  uniquely  determine  the  target  trajectory.  That  is,  the 
functionals  XN{j3)  and  S{0)  are  generally  not  one-to-one.  The  class  of  all  trajectories  jS  such 
that  XsiP)  —  Xn{0t)  or,  equivalently,  E{P)  =  represents  the  inherent  ambiguity 

of  the  radar  problem.  A  problem  of  obvious  importance  is  to  find  outgoing  signals  '(p(t) 
which  minimize  this  ambiguity  class. 

We  have  assumed  above  that  the  return  is  due  to  a  reflection  from  a  single  target.  If 
N  distinct  targets  are  involved,  then  we  can  approximate  the  return  as  a  superposition 

.  (7) 

As  noted,  (7)  is  an  approximation  because  it  ignores  multiple  reflections.  Although  these 
can  often  be  ignored,  they  can  also  cause  resonances  (ringing),  hence  must  sometimes  be 
taken  into  account.  This  can  be  easily  done,  in  principle.  For  example,  the  signal  received 
by  the  doubly-reflecting  path  a  —>■  /3m  —*■  /3n  7  Is 


which  can  be  immediately  converted  to  a  triple  integral  by  using  the  definitions  (2)  and 
(3).  Sums  of  contributions  from  various  “trial”  scattering  paths  may  then  be  matched  with 
the  actual  return,  defining  a  generalized  ambiguity  fimctional 


X{Pi,/32,---,Pn)  =  (X,#i./32,..../3n  )> 

and  the  Schwarz  inequality  may  be  used  as  in  the  case  of  a  single  path  to  optimize  the 
match. 

This  method  is  reminiscent  of  Feynman  diagrams^ ,  where  fundamental  processes  are 
represented  by  multiple  integrals  with  corresponding  intuitive  diagrams.  Because  the 
physics  is  built  into  the  formalism  from  the  beginning  through  the  Green  functions,  our 
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model  can  handle  such  complications  in  a  conceptually  straightforward  (if  computation¬ 
ally  nontrivial)  way.  The  resemblance  to  Feynman  diagrams  is  no  coincidence,  and  the 
present  formalism  may  be  modified  to  include  quantum  (photonic)  aspects  of  radar  simply 
by  using  Feynman  propagators  in  place  of  the  retarded  Green  functions. 
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1.  Generalized  Radar  Analysis  with  Point  Wavelets 

2.  Reduction  to  Time-Scale  Analysis 

3.  Reduction  to  Time- Frequency  Analysis 

4.  Directivity:  Radar  with  Pulsed-Beam  Wavelets 

No  previous  knowledge  of  wavelet  analysis  or  time- 
frequency  analysis  will  be  assumed.  Both  will  be  seen 
to  emerge  naturally  from  generalized  radar  analysis! 

COURSE  SUMMARY 

Problem:  Estimate  arbitrary  target  motions  given 
arbitrary,  independent  motions  of  transmitting  and 
receiving  antennas.  Method:  Match  measured  re¬ 
turn  with  trial  return  based  on  trial  target  motions. 
This  gives  the  ambiguity  functional,  maximized  when 
the  two  returns  coincide.  For  monostatic  radar  with 
uniform  target  motion,  it  reduces  to  the  wideband 
ambiguity  function  (time-scale  transform)  of  the  re¬ 
turn.  For  narrowband  signals,  it  further  reduces  to 
the  usual  time-frequency  transform.  To  obtain  direc¬ 
tivity,  use  pulsed  beams.  Time-frequency  anal¬ 
ysis  is  an  approximation  to  wavelet  analy¬ 
sis,  which  is  an  approximation  to  “physical” 
wavelet  analysis. 


LECTURE  1 

Physical  Wavelets  and  Ambiguity  Functionals 

1.  Emission,  Reception  and  Reflection 

For  simplicity,  we  ignore  polarization.  Since  our  mo¬ 
tivation  is  radar  and  remote  sensing,  we  use  the  lan¬ 
guage  of  “antennas”  even  though  we  are  really  doing 
acoustics  here.  The  concepts  easily  extend  to  EM. 

Space-time  notation:  (x,  t)  =  x  G  R'^ 

Time-dependent  source  distribution:  J(x,  t)  =  J(x) 

The  resulting  space-time  wave  F{x)  =  F(x,  t)  satis¬ 
fies  the  wave  equation 

-  Vl)  F(x)  =  J(x) 

in  free  space.  The  causal  solution  is  given  by 


F(x)=  /  G(x-x')J(x')dV, 

Jk* 

where  G  is  the  retarded  Green  function 
G{x)  —  G(x,  i)  A^\^\ 


Now  consider  a  point  source  moving  along  an  arbi¬ 
trary  trajectory 

X  =  r(r). 

For  simplicity  denote  the  corresponding  trajectory  in 
space- time  by  a: 

a{T)  =  (r(r),r)  G  Rl 

While  moving,  the  source  emits  a  time  signal  ‘0(t). 
This  moving  point  transmitter  is  (h^scribed  by  a 
source  function  concentrated  on  a(r  ),  nanu^ly 

/oo 

dr  b^'^^x  —  a{T))v(T). 

-OO 

The  emitted  wave  is  then 


/OO 

dr  G(x  —  a(T))f( 7) 

-OO 


dr  G(x  —  r(r).  t  -  r]i.-(T). 
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For  example,  if  r{t)  =  0,  then  'l’a{x)  is  the  spherical 
wave 


|x|/c) 

47r|x| 


^(i  -  r/c) 
inr 


as  shown  below. 


t 


A  physical  wavelet  with  stationary  source 


If  r(T)  =  ro-f  VT,  then  a  translated  (delayed)  . 

and  Doppler  scaled  version  of  xp  [1].  However,  the 
amount  of  scaling  depends  on  the  direction  of  the 
observer  relative  to  the  velocity  v  of  the  source,  as  is 
intuitively  and  geometrically  clear. 


Translations  and  scalings  are  the  “raw  material”  of 
wavelets,  hence  for  general  a(r)  we  call  ^a(^)  the 
acoustic  wavelet  generated  by  xP{t)  along  a{r)  [1- 
3], 

EM  wavelets  are  defined  similarly  as  waves  generated 
by  moving  elementary  currents. 

Here  “elementary”  means  point  sources,  but  this  can 
be  generalized  to  extended  sources  [4]. 


Given  a  trajectory  a,  define  the  emission  operator 
by  {EaTp){x)  =  i.e., 


Ea  converts  any  given  time  signal  to  the  emitted 
space-time  wave  '^a{x)  {one  to  four  variables). 
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Reception  works  the  other  way:  A  receiver  moving 
along  a(r)  converts  a  space-time  wave  F(x)  into  a 
time  signal  by  measuring  the  wave  along  a(r).  This 
defines  the  reception  operator 


mitter: 

Fiefi{x)  =  {Ei3R0F){x) 

/OO 

dr  G{x  -  I3{t))F{P{t)). 

‘OO 


(3) 


{RaF){T)  =  go^Fiair)) 


(2) 


where  is  a  gain  factor.  Thus  Ra  converts  functions 
of  four  variables  to  functions  of  one  variable. 

More  general  reception  operators  are  possible.  For 
example,  in  EM,  we  can  measure  the  induced  current, 
which  involves  derivatives  of  the  field.  The  above 
model  is  the  simplest. 


This  is  similar  to  the  EM  case,  where  the  incident 
wave  induces  a  current  on  the  target,  which  in  turn 
radiates  the  scattered  wave. 

Problem:  Given  a  transmitter  following  an 
arbitrary  trajectory  a(r)  and  a  receiver  fol¬ 
lowing  an  arbitrary  trajectory  j(t).  Given  the 
transmitted  signal  iP{t)  and  the  received  sig¬ 
nal  Assuming  there  is  a  single  target, 

estimate  its  trajectory  Pt{x). 


2.  The  Ambiguity  Functional 


Now  consider  a  point  target  moving  along  an  arbi¬ 
trary  trajectory  /3{t)  through  a  wave  F{x). 

\\e  model  the  reflected  wave  Frefi{x)  by  assuming 
that  the  target  acts  first  as  a  receiver,  then  as  a  trans- 


Let  j3{T)  be  a  trial  target  trajectory.  Then  emission 
of  xp  from  a,  reflection  from  P,  and  reception  at  7  are 
given  by 

xp  —*  EaXp  — >  EpRpEaXp  — >  RyEffRffEaXp. 
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Hence  the  trial  return  at  7  is 

ijjp  —  RyEffRpEoi'ip^ 

where  only  the  dependence  on  the  trial  trajectory  0 
has  been  made  explicit.  Of  course,  ipp  also  depends 
implicitly  on  the  icnown  trajectories  a  and  7. 

Inserting  the  definitions  (1)  and  (2)  gives 

=  9y90  JJ  C'(7(r)  -  0{t')) 

X  G{0{t')  -  a{T"))il}{T"). 

This  is  the  computed  trial  return  eissuming  a  trial 
trajectory  0. 


^'is  emitted  by  a,  reflected  by  0,  received  by  7. 
To  estimate  0,  match  xpp  with  the  actual  return  x: 


{X,'>P0)  =  J 


dr  x{r)  tppir), 


where  we  assumed  that  ipir)  is  real,  hence  also  ipp{T) 
and  x(''')-  The  right-hand  side  depends  on  the  tra¬ 
jectory  0  as  a  whole,  hence  it  is  a  functional  of  0.  We 


call  it  the  ambiguity  functional  and  denote  it  by 
X(/?).  Explicitly, 

X(^)  =  909-y  ///  ^"^  -  Idir')) 

xG{0{T')-a{r"))-il}{T"). 

In  practice,  both  ippir)  and  x{r)  have  finite  energy: 


/OO 

dr  -ippir 

-00 

/OO 

dr  x{r) 

-OO 


r)  <  OO 


Then  x(/3)  =  (XiV'/?)  satisfies  Schwarz’s  inequality: 

\(X,‘>P0}\  <  llxll  \\'<P0\\  for  all  0,  and 
I(X>'0/3)I  =  11x11  llV’/sll  xW  =  Gippir) 

for  some  constant  C. 

Thus  to  estimate  the  actual  target  trajectory  0t  , 
vary  the  trial  trajectory  0  to  maximize  the  normal¬ 
ized  ambiguity  functional 

/m  -  X(0) 


Xn(0)  — 


Then 

|x»m|<  llxll 

|x»(/?)|=  llxll  x(t)  =  CV’u(t). 

We  will  show  that  Xn{0)  is  a  generalization  of  the 
usual  wideband  and  narrowband  radar  ambiguity 
functions. 

Even  better;  Define  the  error  functional 


S(0)  =  1  -  =  1 

Then  Schwarz’s  inequality  reads 


Xn{0) 


m\  ■ 


0  <  <  1  for  all  0,  and 

£{0)=O  x{r)  =  C-ippir). 

To  estimate  the  target  trajectory,  vary  the 
trial  trajectory  0  to  minimize  £{0)-  For  a 
given  0,  we  may  regard  the  value  of  £{0)  as 
an  estimate  of  the  error. 
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Warning:  Even  S{0)  =  0  does  not  guarantee  that 
3  is  the  true  target  trajectory,  only  that  the  trial 
return  is  indistinguishable  from  the  actual  return 
X-  But  clearly  we  cannot  do  any  better  than  that! 
For  example,  in  monostatic  radar  with  point  sources, 
all  target  trajectories  with  the  same  |r(t)|  have  the 
same  S{p). 

It  is  impractical  to  try  all  possible  trajectories  13.  Any 
prior  knowledge  or  reasonable  guess  can  be  used  to 
restrict  the  class.  For  example,  assume  the  target 
moves  with  constant  acceleration  a: 

1  2 

r(T)  =  ro  +  vqt  +  -  ar  . 

Then  £(0)  reduces  to  an  ordinary  function  of  the 
target  parameters: 

£(/3)  5(ro,vo,a), 

which  must  be  minimized  to  get  the  best  parameters 
in  this  class. 


3.  Monostatic  Radar,  Uniform  Target  Motion 

To  see  the  connection  with  ambiguity  functions,  as¬ 
sume 

4  Monostatic  Radar:  The  transmitter  and  receiver 
are  both  at  rest  at  the  origin: 

a{t)  =  7(t)  =  (0,t). 

4k  Uniform,  radial  motion:  The  target  motion  is 
r(t)  =  r(t)n,  where  r(t)  =  |r(t)|  =  ro  +  vt. 
Thus 

/3{t)  =  {{ro  +  vt)n,t) 

is  parameterized  by  two  just  scalars:  The  range  ro  > 
0  and  the  range  rate  v.  If  u  <  0,  the  target  is  closing 
in,  and  if  u  >  0,  it  is  receeding. 
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Theorem:  Under  the  assumptions  4k  and  X,  the 
trial  return  is  given  by 

•ipisit)  =  A{t)ip{at-T),  (4) 

where 

A(t)  =  =  attenuation  factor 

^  ’  167r2r(t)^ 

a  =  - — ^  =  time  scaling  factor  (^) 

c  +  v 

T  =  =  time  shift  . 

c  +  v 


Each  quantity  in  (5)  has  a  simple  interpretation: 

1.  The  attenuation  factor  A{t)  implies  that  the  in¬ 
stantaneous  power  of  Tp/sit)  is 

ipffitf  a  rCt)""*, 


2.  The  scale  factor  a  represents  the  Doppler  effect: 

-c  <  V  <0  =»  a  >  1  and  -tp  is  compressed. 
0<u<c  =>  0<cr<l  and  rp  is  stretched. 


3.  The  shift  r  >  0  represents  the  time  delay  due  to 
the  finite  propagation  speed. 


Proof:  Denote  the  wavelet  received  at  /3(t)  (before 
re-emission)  by  tp'it).  Then 


ip'(t)  =  {R0Ea1p){t) 

=  giij  dt'Gm)-a{t'))^{t') 

=  9^!  dt'  G{r[t)n,t-t')ip{t') 


=  90 


47rr(t) 
ip{t-r{t)/c) 
47rr(t) 


in  agreement  with  the  radar  range  equation. 
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The  reflected  wavelet  received  at  'y{t)  =  (0,  t)  is 
therefore 

V0{t)  =  9i  J  dt'  G(7(i)  - 

=  9tJ  dt' G{-r{t')n,t-t')ip'{t')  (6) 

= 

But 

t-t'  -  =  c(t  -  t')  -  To-  Vt' 

c  c 

ct  —  To  —  {c  +  v)t' 


where 


c  +  v  ct  -  rp 

c  [  c  +  v 

C 


ct  -  To 


6{t  -t'  -  r{t')/c)  =  8  (  ^  {t'  -  T) 


5{t'  -  T). 


Therefore  (6)  gives 


xpffit)  = 


9^90 


C 

C  +  V  47rr*(T) 

c  7/;(r-r(T)/c) 
c  +  ?;  [47rr(r)]2 


r(T)  =  ro  +  u  (  ^ 

\  c  +  v  J  c  +  v  ^  ' 


hence 


j,  HT)  ^(0  _  ct-rp  _  rp  +  vt 

c  c  +  v  c  +  v  c  +  v 

\c+vj  c+v 


Putting  all  this  together  gives 

=  A{t)i}{at  —  t).  QED 


Under  4  and  Jit ,  the  ambiguity  functional  reduces  to 

/OO 

dt  A{t)  -ipiat  -  t)  x{t)  =  x{<^,  t)  . 

-OO 

Below  we  will  see  that  this  is  essentially  the  wavelet 
transform  of  x(i))  also  known  as  the  wideband 
amiguity  function. 


LECTURE  2 

Reduction  to  Time-Scale  Analysis 

Recall  that  for  monostatic  radar  (4)  and  uniform, 
radial  target  motion  (<^),  the  trial  return  is 

lo7r^r{t) 


In  practice, 


r{t)  ^  r(0)  =  ro 


(lurins  the  reflection  time  interval.  Hence 


and  the  trial  return  is 


ip0{t)  =  Ap  -ip{crt  -  t). 
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To  connect  a  with  the  Doppler  effect,  suppose  the 
outgoing  signal  is  time-harmonic  (CW): 

/>0. 

Then  the  trial  return  is 

Mt)  =  ^2.iaft  _ 

The  factor  represents  the  time  delay,  and 

g27rt<7/t  jg  reflected  wave  with  the  Doppler-scaled 
frequency  af.  As  expected, 

—c<v<0  a  >  1  =»  <yf>f 
0<t;<c  =4*  0<o-<l  =>  af  <  f. 


Returning  to  the  general  signal  ip{t),  we  have 


/OO  a2 

dt  V»(crf-r)2  =  — ^  \\ipf. 

-OO  cr 

We  may  assume  that  the  outgoing  signal  4^{t)  has 
unit  energy,  i.e., 


/OO  a2 

dt  i>{tf  =  1,  =>  ^  . 

■OO  0- 

Then  the  normalized  ambiguity  functional  reduces  to 

IW«II 

/OO 

dt  'tp{at  —  T)x{t) 

-OO 

=  (  '0(7, T  )  X  )  ) 


Only  in  the  narrowband  approximation  can  f  ^  af 
be  replaced  by  a  uniform  If equency  shift  f  ->■  f  +  <p- 

where 

7p^,r(t)  =  y/aip(at-T) 

is  a  scaled  and  translated  version  of  ip{t)  with  unit 
energy.  Note  that  H’l/'a.TlI  =  Ill'll  =  1- 


The  set  of  functions 

{■^CT.T  :  cr  >  0  and  r  =  any  real  number  } 

is  called  the  wavelet  family  generated  by  'tp,  and  tp 
is  called  the  basic  (or  “mother”)  wavelet. 

The  normalized  ambiguity  functional  thus  reduces  to 
an  ordinary  function  of  the  target  parameters  (cr,  r) , 


Our  assumptions  and  X  say  nothing  about  the 
outgoing  signal  Tpit).  In  particular.  ip{t)  can  be  an 
ultra-short  pulse  or  an  ultra-wideband  signal. 

Thus  in  the  radar  literature,  x(o‘,  t)  is  known  as  the 
wideband  ainbiguity  function  of  the  return  x(t)- 

In  fact,  the  parameter  cr  has  a  meaning  even  when 
ip{t)  is  an  ultra-short  pulse,  since  the  "backscattering 
transformation” 


Xn{0) 


x(o->  t)  =  { 


.x) 

OO 

dt  'tp{at  -  T)x(t) 

*oo 


which  is  called  the  wavelet  transform  of  x(t)  with 
respect  to  the  wavelet  family 

Note  that  for  radar,  only  with 


2ro  ^  1  c-v 

T  = -  >  0  and  cr  = -  «  1 

c+v  C+V 

are  needed.  This  will  also  be  understood  in  the  con¬ 
text  of  wavelet  analysis.  For  now  let  r  be  any  real 
number  and  cr  >  0. 


'tp{t)  -*  tpaA^)  =  Vaipif^t  -  T) 

occurs  entirely  in  the  time  domain. 

When  ■0(f)  is  narrowband,  we  will  s<‘(‘  that  the 
Doppler  scale  factor  a  can  be  replaces  1  by  a  Doppler 
frequency  shift.  For  this  reason  we  call  (cr.  r)  the 
wideband  target  parameters. 

The  general  method  of  ambiguity  func  tionals  and  er¬ 
ror  functionals  now  takes  the  following  reduced  form. 

The  Schwarz  inequality  gives  ’ 

|x(cr,r)|=|(0a,r,x)|<  ll'0<T.rll  !l\il  =  !kll, 
|x(<2-,t)|=  IIxII  x{i)  =  Cvo.rit)- 
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To  find  the  best  wideband  parameters,  we  must  there¬ 
fore  maximize  |x(cr,r)|. 

Under  and  A,  our  physical  wavelet  analysis  of  the 
return  by  its  ambiguity  functional  reduces  to  ordi¬ 
nary  (“mathematical”)  wavelet  analysis  by  the  time- 
scale  wavelet  transform! 

Altenatively,  we  can  define  the  wideband  error 
function 

^)  =  i  -  T-W 

ll^..rP||xiP  lIxlP  ■ 

Then  the  objective  is  to  minimize  £{a,T),  since 
0  <  £{(j,t)  <  1 

£((7,r)  =  0  <;==>  x{i)  =  Cipa,Tit)- 

Furthermore,  the  value  S((7,t)  gives  an  estimate  of 
the  error  in  the  trial  return 

Once  the  wideband  parameters  (cr,  r)  are  estimated, 
the  range  and  range  rate  are  given  by 


1  +  cr 


and  V  = 


To  understand  the  connection  between  time-scale 
analysis  and  time-frequency  analysis,  look  at  the 
wavelet  transform  in  the  frequency  domain. 

We  use  the  Fourier  transform  notation 


X(/)  = 


3-27r  ift 


Xit)=  /  df  e2-/‘x(/). 


/oo 

dt  s/ail){at  -  r) 

-CXD 


g— 27r  ifrfa 


fpif/cr)- 


Parseval’s  identity  states  that  the  inner  product  in 
the  time  domain  equals  the  inner  product  in  the  fre¬ 
quency  domain,  hence  (since  is  complex) 

X(£T,  r)  =  (  V’tr.r  ,  X  )  =  (  ^<T,T  ,  X ) 

/OO 

df 

-OO 


Inserting  (1)  gives  the  frequency  domain  version 
of  the  wavelet  transform: 


9 


We  require  ip  to  satisfy  the  oscillation  condition 


(O) 


That  is,  At)  has  no  “DC  component”. 

This  is  physically  reasonable:  Since  the  DC  compo¬ 
nent  is  constant  (static),  it  cannot  propagate  and  may 
as  well  be  removed.  In  wavelet  analysis,  a  require¬ 
ment  closely  related  to  (O)  is  called  the  admissibil¬ 
ity  condition. 


If  ip(t)  is  integrable,  then  ip(f)  is  continuous  and 
ip{f)-*0  as  f-^0. 

Since  ip{t)  has  finite  energy,  we  also  have 
ip{f)  -*0  as  I/I  ^  00. 

Suppose,  in  fact,  that  Af)  is  a  bandpass  filter,  i.e., 
ipif)  ~  0  outside  the  double  band  /i  <  |/|  <  /2 
for  some  given  frequencies  0  <  /i  <  f 2-  Then 

-7=  Af  A)  ~  0  outside  of  afi  <\f\  <  cr/2  . 
yja 

That  is,  a~^l'^Af  A)  is  a  scaled  bandpass  filter. 
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We  can  thus  interpret  ^  as  follows: 

1.  Localize  the  spectrum  x{f)  using  the  scaled  filter: 

Xif)  -7=  Xif)  =  X<r(/)- 

n/ct 

2.  Take  the  inverse  Fourier  transform  to  get  x{(^,t): 

rOO 

xAf)^  df 


The  wavelet  transform  x  can  be  inter¬ 
preted  35  a  scale-windowed  inverse 
Fourier  transform  (SWIFT)  of  x- 


In  the  narrowband  approximation,  this  view  will  con¬ 
nect  the  time-scale  analysis  smoothly  with  the  usual 
time-frequency  analysis,  where  the  window  slides  in¬ 
stead  of  being  scaled;  see  also  [3] . 


But  what  can  wavelet  analysis  contribute  to  radar? 
The  wavelet  inversion  formula,  for  one. 


•  Given  only  x{cr,T),  find  xi^)- 


This  can  serve  to  model  a  dense  target  environment. 
To  find  x{t)  from  x(o-,t),  note  that  says 

x(a,T)  ’^(//<7)*X(/)|  (r/cr). 

Taking  the  inverse  Fourier  transform  with  respect  to 
r/cr  gives 

/oo 

drer2-^"/-X(<T,r). 

-OO 


Multiply  both  sides  by  cr  jc)  and  write 

for  the  spectrai  energy  density  of  V’-  Then 


5,/,(//(t)  xif) 


yja 

dr  ^<7,t(/)x(<^,  t). 
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Now  integrate  over  cr  with  an  arbitrary  weight  func¬ 
tion  u;(cr)/cr: 

J  ^  w{a)S^{f/a)xif) 

=  r  drrP,Af)xicr,r). 

Jo  ^  J  — oo 


Define 


Then  w(t)  is  the  impulse  response  of  W  and 

/oo 

du  w{t  -  u)  x(u) 

-oo 

/oo 

df  e2-/‘W(/)x(/). 

-OO 

Taking  the  inverse  Fourier  transform  of  (3)  gives 


rOO 

W{f)  =  /  —  wia)S^{f/a),  (2) 

Jo 

the  frequency  filter  corresponding  to  lo(cr).  Then 

l'F(/)x(/)  ~  [[  dT  ‘wio’) 'Ipcr^rif)  Xi^J '^)' 

(3) 

The  inverse  Fourier  transform  of  the  left  side  is 
the  convolution  operator  W  (“stationary  system”) 
with  transfer  function  W{f).  That  is,  let 

/OO 

df  e'^^^^^Wif). 

-OO 


W  * 


x(t)=  r  df  e2-^‘W(/)x(/) 

J— OO 

^IL 


—  dr  w{a)  i>c,Tit)  x(o")  t)- 

<r>0  ^ 


0 


The  meaning  of  0  is  as  follows: 

We  interpret  w(t)  as  the  impulse  response  of  a  filter. 

Then  C>  shows  that  the  filtering  can  be  done  in  the 
scale  domain  instead  of  the  frequency  domain,  with 
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the  correspondence 


This  is  the  reconstruction  formula  for  the  WT. 
A  necessary  condition  is 


I  frequency  filter  w{a)  =  scale  Biter 

j  Fourier  basis  =  <->•  =  wavelet  basis 

\  FT=x{f)--xi<T,r)  =  WT. 

O  gives  the  reconstruction  of  x{t)  from  x{<r,  r)  as  a 
special  case:  Choosing  w{a)  =  1  gives  (for  /  7^  0) 


w{f)  =  r  ^  smm = /°°  f  SAC 

Jo  <y  Jo  ^ 

“i”f 

a  constant.  Hence  w  *  x(^)  —  0  gives 
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the  computation.  Any  other  a  priori  knowledge 
about  the  target  motion  can  be  used  similarly 
to  “compress”  (fZ). 

But  note  that  in  sonar,  the  relative  velocity  of  the 
target  need  not  be  negligeable  compared  to  the  prop¬ 
agation  speed  of  the  acoustic  waves.  (Consider  a 
torpedo  and  a  ship  closing  in  at  maximum  speeds.) 
Then  the  reconstruction  formula  (1Z)  still  works,  with 
the  proviso  that  larger  scale  factors  become  relevant. 
By  contrast,  the  narrowband  ambiguity  function  (dis¬ 
cussed  in  Section  2  of  this  lecture)  fails  to  take  such 
situations  into  account. 

Applications  of  (7Z)  and  0: 

1.  Dense  Target  Environment:  A  single  tar¬ 
get  with  wideband  parameters  (cr,T)  gives  a  return 
tiV.r(0- 

Therefore  we  model  a  dense  distribution  of  targets 
(e.g.,  clouds  or  dust  in  wind)  by  a  “target  density 
function”  T)(cr, r),  with  return 

Xit)  =  ff  —  dripa^rit)  D{a,T). 

J J a>0  ^ 


=  Jo  Y  (■^) 

called  the  admissibility  condition  for  ip. 


Note  that  A  =7>  ip(0)  =  0,  the  oscillation  condition! 


In  (7Z),  all  delays  t  €  R  and  all  scales  cr  >  0  are 
used.  But  on  physical  grounds,  we  expect  only  r  >  0 
and  a  w  1  to  be  relevant,  since  \v\  c  in  radar. 
This  is  reflected  in  (7^)  by  the  fact  that  for  negative 
T  or  large  |logcr|,  the  coefficient  function  x((T,  r)  is 
negligeably  small  since  ipa.rit)  is  a  poor  match  for 
the  actual  return  x(^)- 

In  fact,  we  can  restrict  the  integral  in  (TZ)  to  the 
values  (cr,  r)  expected  to  be  most  likely,  simplifying 


The  objective  is  to  solve  for  D(<t,t)  in  terms  of  x- 
[fZ)  suggests  that 

D(o,r)  =C^1x(<^,t), 

but  this  is  not  quite  right.  A  deeper  analysis  shows 
that  D{cr^T)  cannot  be  found  using  a  single  wavelet 
'ip.  Instead,  a  whole  class  of  different  basic  wavelets 
(each  with  their  family)  must  be  used  to  estimate 
D(cr, r).  See  [3,  5]  and  references  therein. 

2.  Doppler  Scale  Processing;  We  can  process 
the  WT  of  the  return,  by  choosing  a  weight 

function  w{a)  and  using  <).  This  might  be  done  to 
remove  noise  or  to  concentrate  on  a  particular  range 
of  scales,  refining  the  constraint  a  ^  1. 

Since  a  is  directly  related  to  velocity,  such  scale  fil¬ 
tering  can  be  used  for  moving  target  ID  and  clut¬ 
ter  rejection.  In  the  wideband  case,  this  should 
work  better  than  Doppler  frequency  filtering,  since 
the  Doppler  frequency  is  frequency-dependent! 
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LECTURE  3 

Reduction  to  Time- Frequency  Analysis 
Sampling  in  Time-Scale  versus  Time-Frequency 


The  resulting  time  signal 


1.  The  Narrowband  Approximation 

Under  what  conditions  can  Doppler  scaling  be  re¬ 
placed  by  a  Doppler  frequency  shift?  Again  let  ijj  be 
a  bandpass  filter: 

'ip{f)  ~  0  unless  0  <  /i  <  |/|  <  A 

for  some  0  <  /i  <  /2  <  oo.  Then  positive  and  nega¬ 
tive  frequencies  move  in  opposite  directions.  For  ex¬ 
ample,  suppose  that  a  >  1  (target  closing  in).  Then 


/>0 

/<0 


<yf  >  f 
<  /• 


Thus,  in  order  to  approximate  scaling  by  shifts,  we 
must  eliminate  the  negative  frequencies  in  ^(/): 


^(/)  ^'^(/) 


^U(/)  if/>0 
“  I  0  if  /  <  0. 
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is  complex.  It  is  called  the  (Gabor)  analytic  signal 
of  4)  because  it  extends  analytically  to  the  upper-half 
time  plane.  (See  [3],  Section  9.3  for  a  generalization 
to  many  variables,  like  space-time.) 

All  the  frequency  components  in  move  in  the 

same  direction  under  scaling  (to  the  right  if  cr  >  1,  to 
the  left  if  0  <  cr  <  1.)  But  their  motion  still  does  not 
look  like  a  uniform  shift,  unless  the  endpoints  move 
by  approximately  the  same  amounts.  To  ensure  this, 
we  must  assume  that  satisfies  the  following 

condition: 

The  frequency -band  [/i,/2]  is  narrow,  in  the  sense 
that 

(AT) 
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Then  we  can  define  the  carrier  frequency  and  the 
bandwidth  of  V’  as 

j,  /z  +  /i  j  I  _  /z  “  /i  ^  r 

fc  =  —7^ —  and  h  =  — ^ —  <  /c  • 

can  now  be  regarded  as  a  slowly  varying  “mes¬ 
sage”  imprinted  on  a  rapidly  oscillating  carrier  wave. 
In  practice,  oscillates  much  too  rapidly  (at  RF) 
to  be  sampled,  hence  the  carrier  is  removed  by  de¬ 
modulation: 

^(t)  =  e-2"»^<=*t/;+(t). 

This  is  the  complex  envelope  or  the  video  signal  of 

Ip. 

Its  real  and  imaginary  parts  are  called  the  in-phase 
and  quadrature  components: 

I(t)  =  Re^'(f)  =  in-phase  component 
Q{t)  =  Im^(t)  =  quadrature  component  . 

A  given  frequency  /  in  [/i,/2]  is  shifted  by 

2u  ,  2u  , 

Af  =  af-f  =  {(r-l)f  = - —  /  « - /, 

■'  J  J  \  '■>  c-\-V  c 


since  |7;|  c.  If  the  narrowband  aissumption  (A/") 
holds,  we  can  substitute  fc  for  /  on  the  right  to  obtain 


where  Ac  =  c/fc  is  the  carrier  wavelength. 
The  point  is  that  now  the  frequency  shift 


(1) 


is  independent  of  f. 


In  the  narrowband  approximation,  all  frequency 
components  in  tp'^if)  undego  a  uniform  Doppler 


shift  (p  given  by  (1). 


The  wideband  parameters  {a,  r)  are  thus  replaced  by 
the  narrowband  parameters  These  are  re¬ 

lated  to  the  target  parameters  {ro,v)  by 


and  r 


3-3 


3-4 


where  we  have  used  |u|  -C  c  in  the  second  equation. 
The  new  objective  is  to  estimate  0  and  r,  from  which 
we  obtain  the  range  and  range  rate  by 

CT  ,  Ac0 

ro  =  y  and  v  =  — —  . 

It  can  be  shown  [2]  that  in  the  narrowband  approx¬ 
imation,  the  wavelets  based  on  the  analytic  signal 
become 


where 


^ At), 


^4, At)  =  -  r)  (2) 


is  the  video  signal  of  xjja.r-  If  is  a  translated  and 
modulated  version  of  '^'{t).  .The  scaling  has  now  been 
replaced  by  a  modulation. 

Equation  (2)  states  that  'if^,T{t)  is  a  phasor  rotating 
at  the  Doppler  frequency  <p.  This  explains  why  the 
complex  video  signals  must  be  used.  If  only  the  in- 
phase  or  the  quadrature  components  were  used,  then 


the  phasor  would  be  replaced  by  a  harmonic  oscillator 
moving  in  a  line,  like  — ►  cos(27r0().  Although 

this  is  enough  to  determine  \4>\  and  hence  luj,  it  does 
not  determine  the  sign  of  u,  hence  the  direction  of  the 
target  motion  (towards  or  away  from  the  radar  site). 
Clearly  this  is  a  very  valuable  piece  of  information! 

Something  similar  to  the  mathematical  transforma- 


At) -*  i’'*' {t) {t)  =  e  •0+(t) 

must  actually  be  performed  at  the  hardware  level  in 
order  to  “slow  down”  the  signals  so  they  can  be  sam¬ 
pled  and  processed.  See  Stimson  [7]  for  a  detailed 
explanation  with  wonderful  illustrations. 

The  video  signal  of  the  measured  return  x{t)  is 


Xit)  =  e 


xAt), 


where  is  the  analytic  signal  of  x(t). 


Claim:  In  the  narrowband  approximation,  the  wide¬ 
band  ambiguity  function  becomes  the  windowed 
Fourier  transform  of  X{t)  with  respect  to  'lt(t)  as 
the  basic  window. 

Recall  that  the  wideband  ambiguity  function,  ex¬ 
pressed  in  the  frequency  domain,  is 

X((7,  r)  =  (  A,t  ,  X )  =  ( A,r  ,  X ) 

1 

=  ^  df  e2-/"/‘^0(//(T)*x(/) 


where 


-j 

r{a,T)  =  4=  /  df 

Jo 

r(a,r)  =  4=  f  df 

0-00 


=  v/?/  df 

Jo 

roo 

=  v^/  df 

Jo 


A-fAYx{-f) 

Af/<^)AfT 


x((t,t)  =  2Re  x‘^((T,r). 

But  x''’(^i'^)  be  expressed  directly  in  terms  of 
the  analytic  signals  and  x'*'  by 

X*A,r)  =  {A,T,X^) 

=  {A,r,xA-^{^o.r.X). 

where  the  last  relation  holds  in  the  narrowl)and  ap¬ 
proximation.  Note  that  the  carrier  !)ha.'<()r 
which  is  contained  in  both  A  AO  \  *(0.  is  auto¬ 

matically  demodulated  by  matching  the  two  analytic 
signals. 

The  narrowband  ambiguity  function  is  defined 
by  matching  the  video  signal  X{t)  of  the  return  with 
the  video  signal  '^^^r{t)  of  the  trial  return; 

X{<t>,T)  =  {<i!4>.r,X) 


r  dt  -  T)'  X{t). 

J  ^oo 


=  AA,t)*. 


This  is  seen  to  be  the  windowed  Fourier  trans¬ 
form  of  X{t),  using  'I[(t)  as  a  basic  window. 
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Thus  we  have  shown  that  in  the  narrowband  approx¬ 
imation,  the  ambiguity  functions  are  related  by 

X(cr,r)  ^  2Re  X{4),t). 

Since  the  wideband  ambiguity  function  is  a  prototype 
of  time-scale  (wavelet)  analysis,  we  conclude: 

Time-frequency  analysis  may  be  regarded  as  the 
narrowband  limit  of  time-scale  analysis,  where  the 
basic  window  is  the  video  signal  of  the  basic  wavelet. 


X{(P,t)  is  used  to  determine  {(I>,t)  in  the  same  way 
as  x{(y,T)  is  used  to  determine  (cr,r).  Note  first  that 

ll«».r|p  =  ll*f  =  IW*lP  =  5  llrf  =  5- 


The  narrowband  parameters  can  therefore  be  esti¬ 
mated  by  maximizing  |X((/), t)|.  Alternatively,  we 
can  minimize  the  narrowband  error  function 


=  1  - 


ll'&piixip 


since  the  Schwarz  inequality  states  that 


0  <  £{(P,t)  <  1  and 
£(0,r)  =  0  ^  X{t)  = 


Hence  Schwarz’s  inequality  implies 

ix(Ar)i<ii't*,,inixii  =  Timi 

\m,T)\  =  T  ||x||  «=!•  x(!)  =  c^tAt)- 


Summary  of  Ambiguity  Functional  Methods 


The  ambiguity  functional  x(/3)  and  error  functional 
£{(3)  can  be  used  to  estimate  an  arbitrary  point- 
target  trajectory,  with  both  the  transmitting  and  re¬ 
ceiving  platforms  in  arbitrary  and  independent  mo¬ 
tions.  When  the  radar  is  monostatic  and  the  target 
motion  is  uniform  in  the  line  of  sight,  then  xiP) 

£{P)  reduce  to  the  corresponding  wideband  functions 
x(cr,T)  and  £(cr,r),  which  can  be  used  to  estimate 
the  wideband  target  parameters  (cr,r).  If,  more¬ 
over,  the  outgoing  signal  tp(t)  is  narrowband  and  the 
target  speed  lu]  c,  then  the  wideband  functions 
further  reduce  to  the  narrowband  functions  X((/>,  r) 
and  £{cf),T),  which  can  be  used  to  estimate  the  nar¬ 
rowband  parameters  In  this  sense,  physical 

wavelet  analysis  is  a  generalization  of  ordinary  (time- 
scale)  wavelet  analysis,  and  the  latter  is  a  generaliza¬ 
tion  of  time-tequency  analysis.  Radar  provides  the 
glue! 


2.  Time-Scale  vs  Time-Frequency  Sampling 

Actual  computations  are  numerical,  hence  use  sam¬ 
pling.  The  wavelet  transform  naturally  divides  the 
frequency  domain  logarithmcally  into  proportional- 
width  frequency  bands,  as  in  music,  so  that  each 
band  is  sampled  at  a  rate  proportional  to  its  mean 
frequency.  In  the  narrowband  approximation,  the 
widths  of  all  the  frequency  bands  become  equal. 
Hence  in  time-frequency  analysis,  all  bands  are  sam¬ 
pled  at  the  same  rate.  This  multirate  feature  makes 
wavelet  analysis  much  more  efficient  than  time- 
frequency  analysis.  To  show  this,  recall  that  for  real 
^p(t)  and  x(l), 

/OO 

dt  ipicrt  -  T)x(t) 

-OO 

1 

=  4=  /  df 

J-oo 

«/  — OO 
rOO 

=  2Re  /  d/ 

Jo  ■ 
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Suppose,  for  simplicity,  that  tp  is  a  bandpass  filter 
\\dth 

i’if)  ~  0  unless  1  <  |/|  <  2. 

Since  and  x(t)  are  real,  it  suffices  to  consider 
only  positive  frequencies,  i.e.,  look  at  x'^{(7,t).  The 
above  shows  that  depends  on  x(/)  mainly 

in  the  band 

a  <  f  <  2(7. 

Therefore,  to  cover  the  entire  spectrum,  we  sample 
X‘^(<7,  r)  at  the  scales 

am  =  2^,  me  Z  =  {0,±1,±2,±3,---}. 

That  is. 


X'*'  {2^ ,  r)  covers  the  band  Bm  :  2”^  <  /  <  2’^ 


Since  the  width  of  Bm  is 

2771-4-1  _  2^771  _  2^ 

we  sample  x'^{2'^,t)  in  r  at  the  Nyquist  rate  2"*. 

To  recover  all  but  the  DC  component  of  x{t)  from 
the  wavelet  transform  x(c7,r),  sample  x'*'(a,T) 
at  the  discrete  scale-time  points  Om  =  2"^  and 
im,n  =  n/2^,  for  all  integers  m,n. 


Note  that  the  above  sampling  is  minimal,  and  may  in 
fact  be  too  sparce  for  some  wavelets.  See  [3],  Chapter 
6. 

It  does  apply  when  xp  generates  a  wavelet  basis  de¬ 
fined  by 

<f>m,n{t)  =  2^^^xP{2^t  -  n)  =  2^/^xP{2^{t  -  n/2^)). 


As  m  varies  through  the  integers,  all  tequencies  ex¬ 
cept  DC  (f  =  0)  are  covered.  Since  the  DC  compo¬ 
nent  does  not  propagate,  this  is  OK. 


This  includes  all  the  compactly  supported  orthonor¬ 
mal  wavelet  bases  discovered  by  Daubeshies  [8]  and 
their  biorthogonal  variants,  as  described  in  [9]. 


Just  as  a  look  at  the  continuous  wavelet  transform 
indicated  how  sampling  must  be  done  in  the  time- 
scale  plane,  so  will  a  look  at  its  narrowband  limit,  the 
windowed  Fourier  transform,  reveal  how  to  sample  in 
the  time-frequency  plane. 

Recall  that  the  narrowband  ambiguity  function  is 


_  if cT 


/oo 

-CX> 


For  simplicity,  assume  that  the  video  signal  '5  is  con¬ 
centrated  on  the  band  1/|  <  1/2,  i.e., 

4'(/)«0  unless  1/1  <  ^  • 

Then  the  above  shows  that  X{<p,  r)  depends  on  X{f) 
mainly  in  the  frequency  band 

4-1  + 

Thus,  to  cover  the  entire  spectrum,  we  sample 
X{<P,t)  at  the  Doppler  frequencies 

(pm  =  m,  m  E  7,. 


That  is. 


X  {<p,  t)  covers  the  band  B'm  :m  —  ^  </<m-l--. 


As  m  runs  through  the  integers,  the  entire  spec¬ 
trum  X{f)  is  covered.  Since  B'm  has  unit  width,  the 
Nyquist  rate  for  sampling  X  (m,  r)  is  also  unity. 

To  recover  X{t)  from  its  windowed  Fourier 
transform  X((p,T),  sample  X{(p,T)  at  the  dis¬ 
crete  time-frequency  points  (pm  =  =  n, 

for  all  integers  m,n. 

As  in  the  case  of  time-scale  analysis,  this  sampling 
is  minimal  and  may  actually  be  too  sparce  for  some 
windows  See  [3],  Chapter  5. 

Note:  The  width  of  the  scaled  band  Bm  'was  2^,  while 
the  width  of  the  modulated  band  B'm  is  constant. 
This  is  what  makes  wavelet  analysis  so  much  more 
efficient  because  it  allows  high-frequency  bands  to  be 
sampled  at  a  proportionatey  high  rate. 
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Macroscopic  Interpretations  of  (a,  r)  and  (0,r) 


LECTURE  4 


Let  •0(t)  and  ^'(t)  be  concentrated  on  [—T,T]. 


Directivity:  Radar  with  Pulsed-Beam  Wavelets 


Let  be  concentrated  on  [/c  —  b,fc  +  6] 

Let  '!'(/)  be  concentrated  on  [—6,6], 

In  time- scale  analysis,  a  represents  the  passband  of 
’b'a.Af)  and  r  represents  the  pass-time  of 
That  is,  they  represent  the  intervals 

a  < — ^  [<r(/c  -  h),a{fc  +  6)] 

Similarly,  in  time-frequency  analysis,  (p  represents 
the  passband  of  and  r  represents  the  pass¬ 
time  of  That  is, 

< — *■  [(l>  ~  b,(j)  b] 
r^[T-T,T-bT]. 

Thus,  both  sets  (cr,r)  and  (</>,t)  are  actually 

macrosopic  frequency  and  time  parameters. 


1.  Complex-Source  Pulsed  Beams 

Given  a  small  source  centered  at  x  €  R^,  responding 
to  an  impulse  at  time  t.  Let  the  response  at  x'  at 
time  t'  be 


/s:(x',t'|x,t)  =  K{x'\x), 

a  solution  of  the  wave  equation  in  x'  with  a  spatial 
source  distribution  representing  the  source  at  time  t. 

If  the  source  is  a  point  at  x,  then 
•  K{x',  t'  I X,  t)  =  G(x'  —  x,t'  —  t)  =  G{x'  —  x), 


where  G  is  the  retarded  Green  function  for  the  wave 
equation; 


A  simple  model  can  be  given  for  extended  sources 
as  follows: 
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1.  Find  an  analytic  extension  K{x'  \  z)  of  K{x'  \  x) 
complex  space-time  source-locations: 

z  =  {z,u)ec\ 

where 

z  =  X  -I-  iy  €  C^  and  u  =  t  -h  is  6  C. 

2.  Interpret  K{x'  \z)  as  a  function  of  the  real  space- 
time  observation  point  x'  =  {x' ,t')  €  R^. 


Heyman  and  Felsen  [10]  showed  that  when  the 
appropriate  extension  K{x'\z')  is  chosen,  it 
can  be  interpreted  as  a  circular  disk  of  radius 
|y|,  centered  at  x  and  radiating  a  pulsed  beam 
in  the  direction  of  y  at  time  t. 


We  will  see  that  the  imaginary  space-time  vector 
(y,s)  must  satisfy  |y|  <  cs  in  order  for  the  exten¬ 
sion  to  be  defined,  and  fyj  <  cs  in  order  for  it  to 
be  analytic.  Heyman  and  Felsen  usually  choose  the 
minimum  value  of  s,  i.e.,  s  =  |y|/c.  We  will  choose 
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s  >  |y|/c,  so  that  K{x'\z)  is  analytic.  Since  the 
value  of  s  is  related  to  the  duration  r  of  the  pulse  by 
r  =  s  — |y|/c,  leaving  s  free  means  that  all  eight  of  the 
real  variables  x,  y ,  t,  s  have  a  physical  interpretation 
in  terms  of  the  source. 

In  [11]  and  [3],  complementary  extensions  K{z'\x) 
were  derived  and  interpreted:  The  location  x  of  the 
source  remained  real,  while  the  observation  point  be¬ 
came  complex,  i.e.,  x'  ^  z'  =  x'  +  iy' . 

It  was  shown  that  the  eight  real  variables  x' ,y'  €  R"* 
can  be  interpreted  as  phase-space  coordinates,  some¬ 
thing  like  an  eight-dimensional  time-frequency-space- 
wavenumber  analysis  of  the  waves.  (See  Chapter  4  of 
[11],  and  Chapters  9-11  of  [3].) 

The  pulsed  beams  will  be  used  to  generalize  the  am¬ 
biguity  functional  formalism  introduced  in  Lecture  1 , 
which  applied  only  for  point  sources  and  targets,  to 
disks:  The  transmitting  and  receiving  antennas  will 
be  dishes,  and  the  target  will  be  a  circular  disk.  More 
complicated  targets  can  then  be  composed  from  such 
disks. 
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T'o  derive  the  pulsed  beams,  recall  K{x'  \x)  —  G{x^  ■ 
x),  where 


G{x)  =  G(x,  i) 


6{t  -  |xl/c) 


We  will  continue  G{x)  to  G{x  +  iy)  by  extending  the 
numerator  and  the  denominator  separately. 

1.  Extension  of  Ixl: 


X  =  vx  •  X  =  vx 


^  =  a/x2  -  y2  +  2ix  ■  y. 

A  branch  of  the  complex  square  root  must  be  chosen. 

2.  Extension  of  S(t  -  jxl/c):  First  extend  6(t)  to 
6(t  +  is),  then  replace 


S(t  -  |x|/c)  — »  6(t  +  is  —  \/z  ■  z/c). 
The  extension  of  6{t)  is  defined  by 


poo 

6{t  +  is)=  / 

Jo 


if{t-\-is) 


which  is  the  analytic  signal  of  6{t).  For  the  integral 
to  converge,  we  need  s  >  0.  Then 

6{t  +  is)  =  ^  ,  s>0. 

2it  t  +  is 

Note  that  the  boundary  values  of  6  are  related  to  6 
by 

lim^  j^2Re^(t  +  is)  =  6{t). 


G{z)  =  G{z,u)  = 


S(t  +  is  —  jc) 
(47r\/z2) 


87r2(t  +  is  —  vz^/c)  vz^ 


For  6{t  +  is  —  jc)  to  be  defined,  we  need 
Im  (t  +  is  —  ■\/^/c')  >  0,  =>  Im  \/z^  <  cs. 


It  can  be  shown  that  holds  for  all  x  if  and  only  if 

|y|  <  cs. 


This  means  that  the  real  four-dimensional  vector 

y  =  (y.'s)  = 

belongs  to  the  future  cone 

=  {(y,s)  e  R'* :  |y|  <  cs}. 

Thus  G{x  +  iy)  is  analytic  in  the  complex  space-time 
domain 

T+  =  {x  +  iyeC‘^:ye  V^}, 

called  the  future  tube  [3,  11]. 

For  the  above  to  make  sense,  we  must  choose  a  branch 
of  \/^.  In  fact,  the  “branch  points”  are  given  by 

=  x^  -  y^  •+•  2ix  •  y  =  0, 


or 

x^  =  y^  and  x  •  y  =  0. 

For  fixed  y  7^  0,  the  Erst  equation  states  that  x  be¬ 
longs  to  the  sphere  of  radius  |y|  centered  at  the  ori¬ 
gin.  and  the  second  equation  further  restricts  x  to 
the  great  circle  in  the  plane  orthogonal  to  y. 

Given  y  ^  0  and  s  >  |y|/c,  the  extended  Green 
function  G(z)  is  therefore  singular  on  the  branch 
circle 

Cy  =  {x  €  :  |x|  =  ly]  and  x  •  y  =  0}. 

.Nh)w  fix  xo  on  Cy.  Given  any  branch  of  \/^,  follow 
this  branch  while  moving  around  the  circle  of  radius 
F  <  |yl  centered  at  xq  and  lying  in  the  plane  through 
the  origin  orthogonal  to  the  plane  of  Cy.  (The  second 
circle  is  like  a  chain  link  around  Cy!)  Upon  traversing 
this  link,  we  find  that  has  changed  sign,  hence 
we  have  moved  to  the  other  branch! 
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Thus,  to  define  G{z)  uniquely,  we  must  introduce  a 
branch  sheet  Sy  (i.e.,  a  two-dimensional  analog  of  a 
branch  cut),  whose  boundary  is  the  branch  circle: 

dSy  =  Cy. 

G(z)  is  discontinuous  across  Sy.  Different  choices  for 
Sy  give  different  physical  solutions  [10].  Since  we  are 
interested  in  compact  sources,  we  choose  for  Sy  the 
disk 


5y  =  {x  €  :  |x|  <  |y|  and  x-y  =  0}, 

which  will  be  called  the  source  disk.  The  reason  for 
this  name  is  as  follows. 

We  define  the  response  at  the  real  space-time  obser¬ 
vation  point  x'  e  R'‘  to  an  “impulse”  disturbance  at 
the  complex  space-time  point  z  =  x  +  iy  eT+  by 


K{x'  I  z)  =  G{x'  -  z*)  =  G{x'  -x  +  iy). 


(1) 


Note  that  the  right-hand  side  is  well-defined  since 
z^T^  implies  that  x'  —  2:*  €  7!^ . 


K{x'\z)  solves  the  homogeneous  wave  equation  at 
all  x'  =  (x',  t')  €  R"^  except  when  x'  —  x  G  5y  ,  sinoe 
K  is  discontinuous  there.  Recall  that  the  EM  ver¬ 
sion  of  K (x'  I  z)  is  interpreted  as  a  (Coulomb)  poten¬ 
tial,  hence  the  EM  fields  are  given  by  its  derivatives 
with  respect  to  x'.  A  discontinuity  across  Sy  there¬ 
fore  gives  a  delta- function  singularity  on  Sy,  which  is 
interpreted  as  a  source  distribution  concentrated  on 
Sy.  I 

The  two  sides  of  Sy  are  not  equivalent:  The  pulsed  | 

beam  is  radiated  in  the  direction  of  y  and  not  — y.  ■ 

Thus  we  can  interpret  K(x'lz)  as  the  output  of  a 
dish  antenna  centered  at  x  with  orientation  given  by 
y,  emitting  a  pulse  of  duration  (roughly)  given  by  s 
at  time  t.  The  vector  y  therefore  acts  as  a  convenient 
“handle”  by  which  the  orientation  of  the  antenna  is 
controlled. 
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2.  Pulsed-Beam  Ambiguity  Functionals 

We  now  use  pulsed  beams  to  generalize  the  ambiguity 
functional  formalism  from  point  sources  and  point 
targets  to  disk  antennas  and  disk  targets. 

Suppose  that  the’ transmitting  antenna  executes  an 
arbitrary  motion,  including  rotations  and  accelera¬ 
tions.  This  is  parameterized  by 

2  =  a(f}  =  x(t)  -h  iy(t),  where 
x(t)  =  (x(t),t)  e  R'^  and  y(t)  =  (y(t),s)  €  V^. 

Assume  |y(t)|  =  constant  <  cs,  but  the  direction  of 
y(t)  may  vary  to  allow  rotations,  tracking,  scanning, 
etc. 

While  executing  this  motion,  the  antenna  is  fed  an 
input  time  signal  Then  the  output  beam  is 

r  dtK(x'\a(t))m- 

J  —00 

We  call  ^q(x')  the  extended  physical  wavelet 
generated  by  ■0(t)  along  a{t). 


Given  a{t),  define  the  emission  operator  Ea  by 


As  before,  Ea  takes  a  function  of  one  variable  (the  in¬ 
put  signal)  to  a  function  of  four  variables  (the  output 
beam). 

If  the  disk  antenna  is  used  as  a  receiver,  it  converts 
space-time  waves  into  time  signals.  This  defines  a 
reception  operator  Ra.  The  simplest  model  for 
the  received  signal  due  to  an  incident  wave  F(x')  is, 
again, _ 

iRaF){t)  =  gaF{a{t)). 


where  Pa  is  a  gain  factor.  Ra  meiustires  the  field 
along  the  complex  trajectory  a{t).  .Again,  more  com¬ 
plicated  receivers  can  be  formulated  which  measure 
derivatives  of  F  along  Oi{t).  (In  EM.  /?,,  could  mea¬ 
sure  the  induced  current  rather  than  th<'  ti(>ld.) 
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^Since  a(^)  is  complex,  the  value  of  F(x')  at  x'  =  a(t) 
must  be  defined.  We  use  the  analytic-signal  trans¬ 
form  (AST),  which  extends  F  to  complex  space-time 
[3,11]: 


This  is  a  space-time  generalization  of  Gabor’s  ana¬ 
lytic  signals,  reducing  to  the  latter  in  one  dimension 
(time  signals). 

If  F{x')  is  a  solution  of  the  homogeneous  wave  equa¬ 
tion,  then  F[x  -f  iy)  is  analytic  in  the  future  tube  T+ 
[3,11]  . 

We  model  reflection  from  a  disk  as  reception  (by 
the  disk,  acting  as  a  receiver)  followed  by  re-emission: 


Here  Qa  is  re-interpreted  as  a  reflection  coefRcient. 

More  complicated  targets  can  be  patched  together 
from  circular  targets  of  various  radii  and  orientations. 
The  reflection  coefficient  then  becomes  a  function  de¬ 
fined  over  the  target  surface. 

We  now  restate  our  elementary  radar  problem  as  fol¬ 
lows: 

A  disk  transmitter  emits  ^/^(t)  while  following  a{t). 

The  wave  is  reflected  from  a  disk  target  following 
/^t(^)‘ 

The  reflected  wave  is  received  by  a  disk  receiver  fol¬ 
lowing  7(t),  giving  a  return  xi^)- 

All  three  trajectories  a(t}, P(t),j(t)  are  complex, 
hence  they  include  arbitrary  scanning,  tracking,  and 
rotation  as  well  as  accelerations. 


Fren(x')  =  (FaRaF}(x') 


x'la(t))F(a(t}). 


Problem:  Given  the  transmitted  and  received 
signals  '0,x  the  motions  a,  7  of  the  transmit¬ 
ter  and  receiver,  estimate  the  target  motion 
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If  we  knew  we  could  compute  x  from 

Thus,  to  estimate  we  choose  a  trial  motion  0(t) 
and  compute  the  trial  return^  which  we  again  denote 
hy'ipp: 

'ippit)  =  {R^EpRfsEa'ipKt) 

=  9^90 1 f  dt'  dt"  K{^{t)  I  m)  Km')  I  «(*"))  m") 

To  estimate  /?t(0)  match  with  the  actual  re¬ 
turn  x{t)  by  taking  the  inner  product  of  the  two  sig¬ 
nals. 

Again  we  denote  the  result  by  x{P)  and  call  it  the 
ambiguity  functional  of  x  (note  that  now  'tpffit)  is 
complex): 

/OO 

dt  i^0{t)*xit). 

-OO 


The  Schwarz  inequality  implies 


lx(/?)  I  =  l(V’/3,x)l  <  II^/jIIIIxII,  and 
ix(^)l  =  \\i’0\\  lixil  ^  x{t)  = 

Define  the  normalized  ambiguity  functional 


Then 

\Xn{P)  I  <  11x11  for  all  trial  motions  /?,  and 

lX;^(/^)l  ~  llxli  ^  x(^)  ~  C'i/>/3(t). 

Thus,  to  estimate  /3t(0»  must  maximize  Xn{I3)- 

Equivalently,  we  can  minimize  the  error  functional 
defined  by 


€{I3)  =  1  - 


Ix(^)P 


since  the  Schwarz  inequality  states  that 
0  <  S{P)  <  1  and  £{P)  =  0  x(t)  =  C-tpffit). 
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|X,v(/3)|  and  S{^)  attain  their  maximum  and  mini¬ 
mum  values,  respectively,  only  when  the  trial  return 
is  indistinguishable  from  the  actual  return. 

Again,  x(0  =  Cip0{t)  does  not  guarantee  /3(t)  = 
l3T{t),  since  the  return  does  not  uniquely  determine 
the  target  motion. 

Reason:  The  functionals  Xn{^)  one- 

to-one,  i.e.,  X/v(A)  =  Xn{02)  does  not  imply  /?i(t)  = 

m)-  • 

Thus  the  radar  problem  is  inherently  ambiguous. 

The  best  we  can  do  is  find  all  jS  with  |Xn(/?)I  =  llxll- 

This  can  be  called  the  ambiguity  class  of  x(t)- 

Waveform  design  problem:  Find  ^  so  as  to  minimize 
the  ambiguity  classes. 
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The  Schwarz  inequality  may  be  used  again  to  opti¬ 
mize  the  match,  just  as  with  a  single  target. 

This  method  is  reminiscent  of  Feynman  diagrams 
[12],  where  fundamental  processes  are  represented  by 
multiple  integrals  with  corresponding  intuitive  dia¬ 
grams.  Because  the  physics  is  built  into  the  formal¬ 
ism  from  the  beginning  through  the  Green  functions, 
our  model  can  handle  such  complications  in  a  concep¬ 
tually  straightforward  (if  computationally  nontrivial) 
way.  The  resemblance  to  Feynman  diagrams  is  no 
coincidence,  and  the  present  formalism  may  be  mod¬ 
ified  to  include  quantum  (photonic)  aspects  of  radar 
simply  by  using  Feynman  propagators  in  place  of 
the  retarded  Green  functions. 


3.  Multiple  Targets  and  Multiple  Reflections 
We  have  assumed  that  there  is  only  one  target. 

With  N  targets,  we  can  approximate  the  return  as 

4>0i,02,---,0N  ~  “^01  d - 1"  '^0N  ■ 

But  this  ignores  multiple  refections. 

These  can  cause  resonances  (ringing)  and  must  some¬ 
times  be  taken  into  account. 

For  example,  the  signal  received  by  the  doubly- 
reflecting  path  a  — >  dm  —>■  dn  7  is 

Ry  Rp^  Ep^  Rp^  Ed  fp , 

which  is  a  triple  integral 

Summing  over  such  returns  with  k  =  0, 1,2, •••  re¬ 
flections  (where  fc  =  0  is  no  reflection,  applicable  to 
bistatic  radar)  .gives  the  complete  A^-target  trial  re¬ 
turn  TpPi,P2,...,0N  ■ 

To  estimate  di^02,  •  •  •  define  the  A^-target  am¬ 
biguity  functional  by  matching  'tppi,p2,...,pN  with  x: 

x{di,d2,---,dN)  =  {tp0i,02,-,0N  -x)' 
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